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Abstract: We propose a definition of vorticity at inverse temperature /3 for Gibbs states in quantum 
XY spin systems on the lattice by testing exp[—/3i7] on a complete set of observables (“one-point 
functions”). We show in particular that it is independent of the choice of a particular basis. Imposing 
a compression of Pauli matrices at the boudary, which stands for the classical environment, we make 
some numerical simulations on finite lattices, and exhibit usual vortex patterns. 

0. Introduction. 

Consider the quantum XY or Heisenberg spin model for 5 = 1/2 on the 2-D lattice Z2, with 
nearest neighbor interactions. Marmin-Wagner, and Hohenberg theorems tell that Gibbs states, for 
all inverse temperature (5, are invariant under simultaneous rotation of spins (absence of continuous 
symmetry breaking in two dimensions). In the classical case, we know a bit more : although there is a 
unique Gibbs state, with rotational symmetry, which rules out the existence of first order transitions, 
a particular form for phase transition exists, characterized by a change of behavior in the correlation 
functions. For the XY system, it has been described by Berezinskii, and Kosterlitz-Thouless in term 
of topological excitations, called vortices [FrSp]. For the Heisenberg model, we observe higher order 
topological defects, called instantons [BePo]. See also [El-BRol,2] for the classical Kac’s model. 

This paper is a first attempt to answer the natural question : How to define vortices in the 
quantum case ? Consider first a system in finite volume A G Zfi. The Hamiltonians are of the 

form ‘h(A), where $ is an “interaction” between sites in A. For nearest neighbor 

XCA 

interaction, the contributing X are pairs {i,j), and the Hamiltonian of the form 



(id) CA 


possibly adding the linear term ^ cjf. We can also glue to A a finite boundary dA and modify Ha{^) 


to H\\jg\(^) accordingly, accounting for an approximate “external field”. In finite volume A U dA, 
the only (normalized) Gibbs state is given by 


( 0 . 2 ) 




tr(e hHAuaAj\^ 


and called the “canonical Gibbs state”. We shall actually define vorticity at inverse temperature (3 
by decomposing the linear form on a canonical (orthonormal) basis of observables. 
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Indeed, to favour the existence of vortices in finite volume, we have completed A by a hnite 
volume environment DA C where Pauli matrices aj are “compressed” in given directions {6j)j^d\. 
Thus, measurements of the observable “direction of spin” are deterministic on dA, which accounts for 
the ’’classical” aspect of the environment, but the sign of spin still remains a free variable, allowing 
for a certain “chessboard symmetry” of the Hamiltonian, and comparison between ferro- and antifer¬ 
romagnetic behaviors. More precisely, for j G dA we replace aj by aj{6j) = liB.ajHej where Hb. is 
the orthogonal projection in the direction 

One of the main requirements for consistency of our dehnition is to check co-variance of the 
vorticity matrix with respect to the choice of orthonormal basis. 

This paper is organized as follows: 

In Part 1, we dehne vorticity matrices as the decomposition of the Gibbs state in a certain 
orthonormal basis of I-point functions; we call quantum vortices the points of the lattice where the 
vorticity matrix vanishes. 

In Part 2, we study holonomy properties on the Lie group SL(2;R), and show how to define (in 
the ideal continuous limit) the class of holonomy of vorticity matrices through the “non-commutative 
degree”. 

In Part 3, we carry some numerical simulations on constrained quantum anisotropic XY systems, 
and provide pictures of vortices at thermal equilibrium. 

In the Appendix we show that vorticity matrices are dehned intrinsically, i.e. don’t depend, up to 
unitary equivalence, on the choice of a specific orthonormal basis within a simple class of observables. 

To some extend, our approach can be generalized to Heisenberg model, replacing vortices by 
“instantons”, or “skyrmyons” as in Belavin-Polyakov theory. But here the non-commutative calculus 
makes difficult to dehne properly the degree. This will be hopefully investigated in some later work. 

Acknowledgments: We thank A.Messager and S.Shlosman for interesting discussions. This work 
was initiated in Hicham El-Bouanani’s PhD Thesis [El-Bo] at Toulon University, and results in Ap¬ 
pendix were obtained with Renaud Ruamps in his unpublished Master’s Thesis [Ru] at Aix-Marseille 
University, both under supervision of the second author. 

1. Vorticity matrices 

Gibbs state (0.1) for spin 1/2 systems, as a linear form on the C*-algebra of observables 

O = 0jeAudAOj, Oj = A^2x2(C) 

(“quasi-local observables” if we were to consider the thermodynamical limit, ) can be decomposed in 
a canonical basis. The simplest way is to restrict to “one-point functions”, i.e. the set O C O of 
2N X 2N, block-diagonal 2x2 matrices, supported on individual sites of AudA, N = jAudAj. More 
specihcally, let again Dr, C O he a real sub-algebra O, of real dimension AN. 

Example 1: Or is the “canonical” algebra, generated by real matrices (D*)igAuaA, whose all non¬ 
diagonal 2x2 blocks vanish, and all diagonal 2x2 blocks vanish, except this supported on site i that 
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takes values in {(5i, (52, < 54 }, where 


0 0 
0 1 


( 1 . 2 ) 




^2 




^3 





So the family of block-diagonal 2N x 2N matrices with 2x2 entry Sj, 1 < j < 4 at the dth place, 
1 < i < iV 


(1.3) 

gives an orthonormal basis (ONB) of 1-point functions Or,. 

Example 2: Or is the algebra generated by Pauli matrices (0®)igAuaA with diagonal block supported 
on site i that takes values in 

We shall restrict to the canonical algebra, whose generators enjoy the nice property of being real 
matrices. Let also or C o be the algebra of 2 x 2 matrices with real coefficients, endowed with the scalar 
product {A\B) = Tr(i?*yl), which is isometric with R^. By extension, the basis 5 = {(5i, ( 52 , ( 53 , ( 54 } 
of Or will be called an “elementary basis” of Or, since N copies of 5, attached to each site i, give a 
basis ( 0 *)igAuaA,i<j <4 of Or. We say the same thing of any other ONB b = { 61 , 62 , ^ 3 , ^ 4 } of or, 
and of the corresponding basis (R*)ieA,i<j <4 of Or, where R* is defined as in (1.3), with bj instead 
of Sj. Actually, the order of the elements of b matters, so we prefer the matrix notation, namely 

(1.4) ^2^gM4x 4(R) 
is a 2 X 2 block-matrix, where bk is of the form 

(1-5) h = (j3l GAf 2 x 2 (R) 

which we identify with the vector b^ = * ( 6 ^^, 6 ^^, 6 ^^, 6 "^^). Actually we will never use the algebraic 
structure of or in this paper. 

Consider now 6 as a linear operator on /Ci 0 JC 2 , with /Ci = /C 2 = R^, and recall [Si,Sect.II.l] the 
partial trace tri from £(/Ci ( 8 ) /C 2 ) to C{JCi) is defined by the requirement : 


trK:i(A(tri( 6 ))) = tr;Ci(g)/C 2 ® 1)6), X £ C{ICi) 


where tr = ^ Tr, d is the dimension and Tr the ordinary trace. Then with the notations above, we 
see easily that : 


( 1 . 6 ) 



tTb2\ 

tr 64 J 


So tri are the components, in some matrix representation, of the usual trace, or (“tracial state”) on 
Or. For simplicity we set T{b) = tri( 6 ) and call it the “matrix of traces”. An important role will be 
played with symmetric basis. 
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Definition 1.1: We call the ONB b symmetric iff the corresponding matrix b in (1-4) is Hermitian, 
i.e. bi = b\, 64 = 64 , and 63 = b^- We call it S-symmetric if moreover b is real, and T{b) has a 
degenerate eigenvalue, that is, is a multiple of identity. We denote by [5]s the class of S-symmetric 
basis. 


Most of the basis are not symmetric, but occasionally we can make them symmetric, by permuting 
or multiplying by —1 some elements. Note that {Id, ia^ ,iay ,ia^} cannot be made symmetric; this 
is one of the reasons why we prefer the canonical basis {5i, <52, (^s, < 54 }. In Appendix, we characterize 
(5-symmetric basis, up to such transformation. 

So far we have constructed “one point functions”, i.e. a basis of Or. In the sequel we content 
with Hamiltonians of type (0.1) which are of second order in the interactions; if we were to include the 

linear term af we would write it as 1^ (8) cj|. Embed Or into Or (8) Or by the usual coproduct 
ieA {i,j) 

A, and set x = A(x) = ^(l0x-l-x0l)eoR0 or, for x G or. So we have “lifted” b = A(b) as a 
family of Or 0 Or by (Bj)i^A,i<j< 4 , with Bj = A(Bj). With the notations of (1.4) and (1.6) we 
have 


( 1 . 8 ) 




V^3 ^4 


G A44Arx4Af(R-), tri(H*) 


f trBl trB2 
V tr B\ tr B\ 


£ A42x2(I1') 


In the same way, we form e Bj, so we can map to each site i G A U dA a 2 x 2 matrix 


(1.9) 




Definition 1.2: 

the matrix : 


We call vorticity matrix at site i, relative to the basis b, at inverse temperature f. 




tri(e ^^H*) 

tr(e-^i?) 


The traceless matrix 


(1.10) H^( 6 ) = H^( 6 )-tr(H^( 6 ))ld 
is called the reduced vorticity matrix at site i. 

Example: A = { 1 , 2 } is a lattice with 2 sites, dA = 0, one has H^((5) = H|(5) = 0. This is observed 
also numerically for all A, with dA = 0, although vortices should merge sponteanously in infinite 
volume. 

If 6 is a symmetric basis of or, then ^^p{b) and Th^j^fb) are hermitean since H is self-adjoint (and 
are real symmetric if H moreover has real coefficients), and 

(1.11) (H^( 6 ))^ = detH^( 6 )Id 

Thus enjoys the nice property, to be diagonalizable with real (opposite) eigenvalues for all sites 

i, and all inverse temperature (3. Viewing these as a field of matrices over the lattice, we can figure 
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out the “vorticity” of the system (whenever this makes sense), by simply looking at their principal 
directions. This also gives a measure of vorticity, i.e. numbers (integers) that should be independent 
of the choice of “elementary” basis b. 

Next we define vortices as the set of sites where the reduced vorticity matrix is singular. 

Definition 1.3: We say that ^ G A is a vortex at inverse temperature /? relative to the 6-symmetric 
ONB b ijf Q^^{b) has a degenerate eigenvalue, i.e. fl|(6) = 0. We call regular the other points. 

By construction, all sites are vortices when /3 = 0. 

Now we turn to consistency of Definitions 1.2 and 1.3 relatively to the choice of b within 5- 
symmetric basis. That 6 is a 5-symmetric basis is a natural requirement for computing the degree of 
D*(6), see Sect.2. With b written as in (1.4), and P G 0(2; R), we set with obvious notations 

(1.12) a = ^PbP 


(i.e. as if fej’s were numbers). The same holds after taking the co-product A of each term, i.e. 
a = bP. This defines conjugacy classes, which pass to the partial traces (1.6), i.e. T(a) = *PT{b)P, 
and T(d) = ^PT{b)P. Moreover, if A G £(R^), we have 


(1.13) 




After lifting a and b to Or (g) Or, (1.8) becomes 


(1.14) 


A* =*P 


Bi 


3 

Bi 


P G M 


4Nx4:N 


(R), tri(A*) = *P 


tv Bl 

tiBi 


tr i ?2 
trBl 


P G A42x2(R) 


and (1.13) also extends when taking X G T(R^'^) and replacing a by A*, b by R*. Let now X = e 
we obtain that conjugacy classes pass to vorticity matrices, i.e. 


(1.15) 


n}{a) = *PD^(6)P, D^(a) = ^P%{b)P 


In Appendix we characterize completely the set of (j-symmetric ONB’s up to the action of the group 
Gq acting on A44(R) by permutations, or multiplication by —1, of some vectors b^ G R^ (recall from 
(1.5) the identification of b^ with the 2x2 matrix bk). Namely, we show in Appendix that if b is S- 
symmetric, then modulo the action of Go, there exists discrete or one-parameter families Pg G 0(2; R) 
such that 


(L16) b = bs = ^Ps6Ps 

So we have : 

Proposition 1.4: Definitions 1.2 and 1.3 are consistent, i.e. vorticity matrices relative to all S- 
symmetric ONB b are related by (1.15) for some Pg G 0(2; R), and in particular is a vortex relative 
to 5 iff this is a vortex relatively to any 6-symmetrie b. 
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Moreover we have the numerical evidence that, among all (5-symmetric basis b, the canonical basis 
6 is most “faithful”, in the sense that have on the boundary lattice dA the same principal direc¬ 

tions as the directions along which Pauli matrices are compressed (associated with the eigenprojector 

n,). 

Let us conclude this section by some heuristic remark: In the case of a translation invariant 
interaction and free boundary condition (5A = 0) it is standard to show that Gibbs state 

(0.1) converges in the thermodynamical limit, and so do the vorticity matrices. This is already of 
interest, because in the framework of Mermin-Wagner and Hohenberg theorems, it is believed that 
(quantum) vortices merge spontaneously in infinite volume, without any boundary conditions. In our 
case however, the compression of Pauli matrices on dA ^ 0 breaks the translation invariance. When 
dA = \ A, we could think of the boundary condition as an external field, that we recall from 

[Si,Sect.II 3] : 

A state is the assignment of an operator px for each finite X C with Tr(/3x) = 1 and 
tr^y(/3xuy) = Px (partial trace on Hy = (8>ieyCf, ) for all disjoint X,Y C T?. Given the nearest 
neighbor interaction ‘h(X), and a state p on Z^, we define the Hamiltonian on all 7 ? 

^aW = - E T^tx\a[{1®Px\a)HX)] 
xnA^0 

that couples A with the external field p through its nearest neighbors at the boundary. 

For any quasi-local observable A on Z^, define the expectation value 

_ TrA(exp[-/?iLP(^)]A) 

^ TrAexp[-/IiL^(ci))] 

Then it is known that both |A|“^ logTrAexp[—/3iL^(<h)] and |A|“^ logTrA(exp[—/3iL^(<h)] A) have a 
limit as |A| —oo, and so has {A)^ On that basis we could expect that the vorticity matrices with 
a conditional external field have a limit in the thermodynamical limit. However, it is not quite clear 
which field p could stand for the compression of Pauli matrices at the boundary (see [AscPil] for the 
1-D XY chain); our Hamiltonian iLAuaA(‘h) is only an approximation for iL^(<h). 

2. Holonomy on the Lie group SL(2;R). 

We pass here to an idealistic continuous limit, where vorticity matrices would be defined as a 
smooth field on (away from vortices) , valued in the Lie algebra sl(2;R), consisting of traceless 
matrices. Our purpose is to integrate such fields vanishing at some points, and define the “non- 
commutative degree”. For advanced results on Differential Calculus on lattices in the scalar case, 
see [Sm]. The non-commutative discrete case, also allowing for an extension of our XY model to 
Heisenberg model, has still to be set up. 

Let M : D C R^ —sl(2;R),x M{x) be a map, such that such that M{x) obeys (1.11), 
i.e. M{xY = A(x)Id, \{x) > 0, and consider p G A^(R^; sl(2; R)) the 1-form defined by 

(2.1) p{x) = ^{M~^{x)dM{x) — dM{x)M~^{x)) 
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with the property of being antisymmetric if M is symmetric. Since 


M ^{x)dM{x) + dM{x)M ^(x)) = 


d\{x) 

\{x) 


we have 


M ^ {x)dM (x) = p{x) + , dM{x)M ^{x) =—p{x) + 


2A(x)’ ' 2 A(x) 

If moreover D C is simply connected and A(x) > 0 in ZD, 2 \(x) exact, and if 7 is a loop in D : 


(2.3) 


M ^{x)dM{x) = — f dM{x)M ^(x) = j p{x) 
J ^ J ^ 


In the general case, using identity M(x)^ = A(x) Id, we find easily 

7 / N dM dM dM dM dM dM dM dxi A dx2 

dp{x) = [M-M- - M-M- - M—M + —M—M] . 

OXi 0x2 0X2 OXi 0X2 OXi OXi 0x2 A^[x) 

and setting M = ^ ^ ^ , a computation shows that 

(2.5) dp = —X~^{adb A dc + bdc Ada + cda A db)M 
Thus the form p is closed if 

(2.6) R{a, b, c) = adb A dc + bdc A da + cda A db = 0 

and this condition holds if M is symmetric. For such a map M uniformly elliptic at infinity, in the 
sense that 

|A(x)| > C > 0, |a:| > ro 


we can define the number 
(2.7) 


Sr^ = det ■ 


27r 


|a:|=r 


p{x) 


which, by Stokes’ formula, turns out to be independent of r > xq. In the same way, if ^ is a vortex 
(i.e. the map M is singular at ^) we define the “local degre local” 

1 


( 2 . 8 ) 


sg = det 


27r 


p{x) 


whenever M(x) is invertible for x 7 ^ integrating on a small contour 7 around 

Let us now compute Maurer-Cartan structure equation for the form p{x) [Ma,p.I65]. The struc- 


1 0 
0 -1 


.62 = 


0 I 
0 0 


.63 = 


0 0 
I 0 


ture coefficients for the Lie algebra sl(2; R), with basis Ci = 

are given hy 2 = “C’l 1 = 2,C'f3 = —= —2,(72 3 = “^32 = 1; C^j = 0 otherwise. For 
i^.pi = E(E Ci jP^ A p’)ek, we find 

k i<j 


— [p,p] = \ ^ {adb A dc + bdc A da + cda A db)M 
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Under (2.6), this relation together with (2.5) show that p verifies 


(2.9) 


dp + [p, p]=0 


Recall that if G is a Lie group, and A its Lie algebra, uj the canonical Maurer-Cartan form on G, 
invariant by left translations, we define Darboux differential of the map / G C^{D;G) by iTf = 
f*uj. The fundamental existence theorem (“Poincare lemma”), with a differential form p G A^(D;^) 
verifying dp + [p, p] = 0, associates (locally) a map / G G^{D-,G), whose Darboux differential is 
precisely equal to p. Moreover this map is unique when assigning its value on a point xq £ D. 

On the other hand we know [Ki,p.ll7&: 321], that the Lie group whose Lie algebra is sl(2;R), 
is the universal covering E = SL(2;R) of the unimodular group SL(2;R). The unimodular group is 
topologically equivalent to the cylinder x R^, its fundamental group equals Z, and E is homeo- 
morphic to R^. (It is known however that one cannot parametrize E by matrices, more precisely E 
cannot be written as a subgroup of some GL(m; C), with m G N, but rather as a tensor product of 
such matrices. ) 

Relation (2.9) ensures the existence of a local primitive N G G^{D; SL(2; R)) of p, the “logarithm” 
of M. If D is simply connected, this primitive is also global. Otherwise, consider its extension to E, 
and let 7 C R^ be a loop at xq G D, we may define the monodromy opeerator acting on functions 
N : j ^ E. 

Example 1: For symmetric matrices in sl(2;R) 


( 2 . 10 ) 


Mq{x) = 


cos nO sin nO 
sin nO — cos nO 


Mr{x) = 


r cos n9 sin nO 
sin n9 —r cos n9 


we have (degf^(M))^ = The 1-form p associated wih Mq is simply 


0 n 

—n 0 


d9. 


Since the fundamental group of Li is Z, s^o and sg are integers so that we set 


( 2 . 11 ) 


Soo = (deg^(M))2 G {0,1,4,9, •••}, sg 


(deg,^(M))2G{0,l,4,9,---} 


Degrees at infinity and at xq are then invariant by homotopy. In particular, a perturbation theory 
can be carried out by expanding M as Fourier series. The degree for matrices M G sl(2; R) verifiant 
(2.6) can be also obtained by Brouwer theory [Mi], by considering these matrices as (locally) a 2-D 
manifold. If matrices M are symmetric, then p is antisymmetric, and 


( 2 . 12 ) 

with n G Z. 



—n 

0 


3. Numerical simulations. 



Recall we have completed the lattice A with an environment dA C where Pauli matrices are 

f cos^ 0 cos 0 sin 0 

compressed in directions i.e. we change a by Healle, where He = 


Thus 


cos y sm 


9 sin^ 9 


(3.1) 


(7f(6»i) = (sin26'i)ne^, (rfifii) = 0 


Hamiltonian (0.1) with nearest neighbor interaction has a large kernel, so it is not directly suitable 
for numerical simulations, even when modified by an external field. In the Classical case, we con¬ 
sider instead the 2-D planar rotator with long range interactions, such as Kac’s model [El-BoRo]; 
renormalizing the Hamiltonian leads to the free energy functional on macroscopic scales (or coarse 
graining), whose critical points are most accessible to numerical analysis. But in the Quantum case, 
renormalization procedures for 2-D planar rotator with long range interactions are not yet available 
(see however e.g. [ScOr] for Ising model). Changing (0.1) to the anisotropic XY model is in fact a 
hrst attempt to lift the degeneracy of the spectrum of the Hamiltonian, and enhance the effects of 
the external field on vorticity even in quite small lattices, so that the predicted vorticity could be 
observed with a fairly good accuracy, only using elementary numerical tools (Wolfram Mathematica 
on a laptop). We discuss below the role of anisotropy. For n,k > 0, consider the Hamiltonian 
(3.2) 




1 

2 (n -|- k) 


{na^ +k af (g) fjJ) 


1 

2 (n -|- k) 


Y, n(af0a^(0,) + a^(0,)0af) 


1 

2 (n -|- k) 


Y nai{9i)(S)a]{9j) 


SO H(^^ j,^{a\dA) is self-adjoint and real. When n = k = 1, A = 7 ?, H = is the most natural 

(isotropic) model with 0"’“(2) symmetry. It enjoys nice properties, like reflection positivity; its spec¬ 
trum is believed to be absolutely continuous on [—2, 2] as this of the Laplacian on Z^, but this is not 
rigorously known, see [DaManTie], [De]. For A: 7 ^ 1, we call / 7 (i^fc)(( 7 | 5 A) the anisotropic XY model. 
Only when dA = 0, is unitarily equivalent to In general, H(^iy.^{a\dA) has no obvious 

symmetry, but it is most suitable for studying vorticity matrices on finite lattices A U dA, at least for 
small /3. 

We consider rectangular lattices of minimal sizes to exclude important volume effects, with suf¬ 
ficiently large dA to constrain the “quantum system” within A. We choose 9j = dojj -|- cj) where ojj is 
the polar angle representing the vector j G dA. We compare calculations for n = 1 and A: = 1, 2,10. 

a) Considerations on speetra. 

Numerically, we observe that the spectrum of H(i^k) is distributed in an interval I close to 
[—2,2], and looks symmetric around 0, allowing for equivalence between ferromagnetic Hamiltonian 
i7(i fc)(cj|9A) and antiferromagnetic —fc)(cj|9A). The distribution has smaller density at the edges 
of I, and larger near A = 0. We present below the integrated density (statistical distribution) of 
states p{X) = #{Afc|Afc < A} in various situations, namely we compare isotropic and anisotropic cases 
without boundary (Fig.l a, b) and anisotropic case with boundary (Fig.l c). 
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tt|A: A<Aol 



Fig.l. The integrated density of states p{\) = #{Afc|Afc < A}: (a) isotropic case without 

boundary, |A| = 19 x 29; (b) anisotropic case i^(i,io) without boundary, |A| = 19 x 29; (c) anisotropic 
case i^(i.io) with 2 boundary layers, |A U 9A| = 23 x 33. 

The reason for degeneracy at A = 0 is the following. In the isotropic case n = k = \ the matrix 
cjf (8 ) ctJ + (8 ) cjJ is of rank 2, so when 5A = 0, half of the eigenvalues of H vanish, and also in the 

general case there is a big degeneracy of the spectrum near A = 0. As in QFT we could try to remove 
that “artificial” part of Ker H by reducing the Hilbert space % = to a “physical space”, but a 
difficulty arises because H is not positive in the form sense. The effect of anisotropy is to lift this 
degeneracy, and enhance vorticity effects. Another reason for degeneracy is boundary effects as can 
be seen from comparison between Fig. l.b and l.c. Degeneracy could be also reduced by enlarging 
the inner lattice, but at the expense of computational difficulties. 

b) Vorticity patterns. 

We study Gibbs state at inverse temperature /3, with significant results provided (3 ranges in 
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some interval, for which however, there is no evidence of a second order phase transition. 

To visualize monodromy of the vorticity matrices n^((5), we plot their principal directions as 
“crosses”, of length proportional to their eigenvalues (recall the reduced vorticity matrices are sym¬ 
metric, with eigenvalues ±^/Xi, ) as we would do with arrows in the classical model [El-BRo]. As 
expected, their principal directions coincide on dA, with those of the eigenprojectors Ilgi.. We expect 
also the number of vortices to be equal to the topological degree d = Soo- Indeed, computing ^ p 
as a discrete integral along a contour 7 G A, not too far from the boundary (in practice, 2 or 3 layers), 
it turns out that the computed degree is close to this we would obtain in example ( 2 . 12 ). 


Because of degeneracy in the isotropic case vorticity matrices in all sites that are not the first 

'c*(/3) 0 


neighbors to the boundary are of the form fl) 3 (( 5 ) 


0 


c*(/3) 


for some c*(/3), and so the reduced 


vorticity matrices are zero fl^((5) = 0 in these sites. This fact can be shown by direct calculations of 
eigenvectors and vorticity matrices and is illustrated in Fig. 2 below. 




I# 





X-V A -V + 

^ X-VA-v 

X'X X X X Jr 


^XXXXX-V-v 

-V X + 


+ A- A- 
+ -V x 

VXXXXXXXXXXXy- + 


(a) 0 


(b) 



Fig 2. Principal directions (lengths correspond to sizes of eigenvalues on a logarithmic scale) for 
different degree: (a) d = 0, (b) d = 1. Lattice is |9A U A| = 15 x 19 with 2 boundary layers. Here is 
considered the isotropic Hamiltonian with /3 = 1. 


If we consider anisotropic case then matrices erf ( 8 ) crj -f erf ( 8 ) erj are no longer degenerate and so 
eigenvalues in the center are no more zeros (see Fig. 3). In this case eigenvalues decay when getting 
far from the boundary towards the center of A and it seems that the rate of decay is exponential 
with distance from the boundary. To make the vorticity patterns more demonstrative we draw the 
crosses on a logarithmic scale. All points sufficiently close to the center of A look like vortices within 
the standard accuracy of computations, but the number of “true” vortices should be equal to the 
topological degree d = s^o- 
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>Ly.xx......*,t-^ 

^*XXXXXXXXX,C^ 




y- ^ 
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Fig 3. Principal directions (lengths correspond to sizes of eigen values in logarithmic scale) for 
different degree: (a) d = Q, {h) d = 1, {c) d = 2, (d) d = 3. Lattice is |9A U A| = 15 x 19 with 2 
boundary layers. Here the anisotropic Hamiltonian lo) is considered with (3 = 1. 


To compute the degree we use the discrete approximation based on finite-differences method: 


(3.4) 


M ^{x)dM{x) 


E 


M{xi 


-1 


M(x 


*+i; 


M{xi 


V'l det M{xi)\J A .^1 detM(xi+i)| x^l detM(xi)| 


and a similar formula for — dM{x)M~^{x). When angles are close to 7 rn/ 2 , n S Z, eigenvalues of 

reduced vorticity matrices are close to zero due to properties of cr* (di) on the boundary. Multiplying 
matrices with a big discrepancy in their eigenvalues would lead to large computational errors; to 
compensate for this effect we use in Eq. (3.4) “normalized” matrices that are divided by square roots 
of their Jacobians. 

The main factor of inaccuracy in degree calculations consists in the discrete approximation of 
the integral and the number of points on an integration contour. As a rule, accurate results require 
a lot of points on the integration contour; but for larger degrees the variation of the angle increases 
from point to point and so do the error due to discrete approximation. This problem can be solved 
by enlarging the lattice size. Fig. 4 presents vorticity patterns with different /3’s for a lattice with 
|5A U A| = 23 X 33 and two boundary layers. 
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Fig 4- The principal directions for different beta: (a) /3 = 1, (b) /3 = 8, (c) ^ = 16 for lattice 
|(9A U A| = 23 X 33 with 2 boundary layers. Lengths of crosses correspond to eigenvalues on the 
logarithmic scale (left) or are set equal (right). Here the anisotropic Hamiltonian iL(i lo) is considered. 
The degree is d = 1. On the figure (a) to the right the eigenvalues in the center that are below the 
computational error are set to zero. 

Another factor that affects the computations is that eigenvalues are decaying while getting further 
from the boundary. This we partially compensate by considering anisotropic model. In principle we 
can consider anisotropic model with smaller anisotropy (e.g. k = 2 oi even close to 1). As soon as 
k > 1 eigenvalues inside the lattice are not zeros and the degree can be computed. But in practice for 
small anisotropy the rate of decay of eigenvalues is very high, so the computed degree is more accurate 
for large anisotropy. To make our vorticity patterns more demonstrative we consider the case of high 
anisotropy with k = 10. In Table 1 below we give the results for degree, computed for a large lattice 
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along the cycle 7 C A consisting of the rectangle of the first or the second neighbors to the boundary 
for different values of the anisotropy parameter = 2; 10. Inverse temperature is /3 = 1. 

Table 1. Table of calculated degree for different values of anisotropy factor k for the first and the 
second neighbors to the boundary. Here the case of 2 boundary layers |Au9A| = 23 x 33 is considered 
with 13 = 1. 


Given k = 2 k = 10 k = 2 k = 10 

degree I'** neighbours 1^* neighbors 2”'^ neighbours 2"''^ neighbours 

1 1.05 1.09 0.89 1.05 

2 1.98 2.03 1.70 1.78 

3 2.76 2.75 2.01 2.50 


Increasing the boundary size for the same inner lattice does not almost change anything inside. 
This can be explained by the “exponential” decay of the information that propagates from boundary 
sites. So the influence of the third boundary layer on the inner points is too small compared with the 
effect of the closer sites on the first and the second boundary layers; in simulations it suffices to use 
only 2 boundary layers. 

Let us consider the influence of the inverse temperature (3. As expected, the smaller j3 the smaller 
are the eigenvalues of when i £ A, because of disorder at high temperature; but taking larger 

(3 makes everywhere on A, very close to which has principal directions (1,1) and 


( 1 ,- 1 ), and the vorticity pattern is destroyed; this is not in contradiction with long range order at low 
temperature, but is probably due to volume effects, (a) 


Let us finally discuss the antiferromagnetic model. It is known that on Z^, the unitary trans¬ 
formation U consisting in flipping the spins at sites i with i odd (i.e. indices i = (^ 1 ,^ 2 ) such that 
Kl = Ih| + K 2 I is odd) intertwines the ferro with the antiferromagnetic models. More precisely, 
—H = U*HU. The reason is that Zg and Z^ (the even and odd lattices) are swapped into each other 
by symmetries on the lines x = n-|-l/2ory = m + l/2 (called the “chessboard symmetry”). There 
follows that tiex.p[/3H]A = trex.p[—/3H]UAU*, and if A = (the canonical basis), we can check 
UAU* = A so the matrices of vorticity (for the Hamiltonian with free boundary conditions) are the 
same. This equivalence holds also in the case of the torus, but not on A C Z^ with an odd number of 
sites. Of course, when dA 7 ^= 0, H and —H are not so simply related; nevertheless, we may observe 
(numerically) that the relation H^(5) = H^^((5) holds with a very good accuracy. 

c) Summary. 

The number of vortices is equal to the topological degree and can be calculated by the integral 
along some contour. We use discrete approximation of this integral to compute the degree for a finite 
lattice A U dA. The main factors responsible for computational errors are the small number of points 
on the integration contour and the fast decay of eigenvalues inside the lattice. We can use larger 
lattices and longer integration contour approaching the boundary to reduce the computational error 
related to the number of points, or consider smaller degrees. We deal with anisotropic case to slow 
down the decay of eigenvalues near the center of A, and take the inverse temperature of order 1 . 
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Appendix 


We provide a proof for the conjugacy relation (1.16) between ^-symmetric basis. Allow for the 
action of Gq, we can assume the matrix 6 E 0(4) with columns (6^, 6^, 6^, 6^) already has determinant 
1, so this is the matrix of a rotation in R^. We recall from [So] the structure of matrices of rotation 
in R^. Matrices of the form 


(a.l) 


AG{a,U) 


( cos a —sin 0*17 

ysinaC/ cosa + sina;j([/) 


with *17 = {x,y,z) a unit vector in R^, and 


/ 0 -z y \ 
j{U) = \ z 0 -X 
\-y X 0 / 


are called left-quaternions. In the same way, 


(a.2) 


A£){a, U) 


f cos a sin a* 17 

y —sinal7 cosa + sinaj(17) 


are called right-quaternions. The point is that any R G 0'’"(R^) can be written as the product of a 
left-quaternion and a right-quaternion. More precisely, there are real numbers a, (3, and unit vectors 
U = (sin 0 cos (^, sin0sin(/), cos 0),V = (sin/ucosz/, sin /r sin z/, cos ^), such that R = Ada, U)Ad{(3, V). 
Thus R is a matrix depending on 6 variables a, /?, 9, (f, /r, u. 

Recall from Definition 1.1 the conditions for 7? to be a 5-symmetric basis, R £ [(5]s; after taking 
some linear combinations between the original equations, this gives the following system: 


(a.3) 

(a.4) 

(a.5) 

(a.6) 

(a.7) 

(a.8) 


sin /3{sin /i[cos a(cos n— sin u) -|- sin a cos 6{cos v -\- sin v)] 

— sin a sin 9 cos ia{cos (j) -\- sin 0)} 

— cos f3 sin a sin 9{cos (/> — sin (/>) = 0 

sin a sin f3 sin 9 sin cos(z/ -|- (/>) = 0 

sin a cos /3 cos 9 -\- cos a sin j3 cos = 0 

sin a sin 0[sin j3 cos cos (j) — cos (3 sin </>] = 0 

sin j3 sin /i[cos a cos v -\- sin a cos 9 sin z/] = 0 

sin a[2 cos (3 cos 0 -|- sin /3 cos{u — (f) sin 9 sin /x] = 0 


This system turns out to be overdetermined, which allows for discrete or 1-parameter families of 
solutions, which can be found with the help of a symbolic computer program. 

Equation (a.4) yields the set of congruence mod tt: (a) a = 0 [tt] or (b) /3 = 0 [tt] or (c) 0 = 0 [vr] 
or (d);U = 0[7r] or (e) zx -|- 0 = |- [tt]. 

In cases (a) or (b), we find a = 0 [vr] and /3 = 0 [vr]. In case (c), a = /3 = ^ [vr] and /x = 0 [vr]. In 
case (d), /a = 0[7r], a = 0 [vr] and /3 = 0 [vr]. So in cases (a) to (d), the solution to the system (a.3)-(a.8) 
consists in a discrete set. 
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Case (e) splits into two sub-cases: (e’) u + (/) = [vr] with a ^ ^ [tt] and /3 7 ^ ^ [tt], or (e”) 
1 / + (j) = ^[tt] with a 7 ^ I [tt] or /3 7 ^ [tt] . 

In case (e’), we find (a = f [vr] and 13 = —^ [vr]), or (a = —^ [vr] and /3 = f [vr]). So in cases (e’), 
the solution to the system (a.3)-(a.8) consists again in a discrete set. 

In case (e”) we find the solutions we obtained before, in addition to a one-parameter set defined 
by a = /3 = 0 = f [vr], 1 / = 0 [vr], and ;U = 0 = t G R. Summing up: 


Lemma a.2: The solutions of R £ [5]^ consist in a discrete set, comprising Rq = ±Id, and 


R^ = i 


/I 

0 

0 




0 

0 


0 

-1 

0 


, Rq = ± 

0 

0 

-1 


0 

0 

-1 

0 

0 

-1 

0 

0 

\0 

0 

0 

1 / 



0 

0 

0/ 


(a.9) 


= zb — 


R^ = zb- 


1 


f 1 1 11 ' 

-1 1 -11 
- 1-1 1 1 
V 1 -1-1 1, 

f 1 1-1 - 1 ' 

- 111-1 
1111 
V-1 1 -1 1 

and 4 one-parameters families, with 0 < s < 1 .• 


Rio{s) = zb 


R 2 = =b — 


, i ?4 = zb- 



Rllis) = zb 

(a.lO) 


Ri 2 {s) = zb 



f " 

^Js{l - s) 

Vs{l - s) 

l-s 

\ 


^ys{l - s) 

—s 

l-s 

-v's(l - s) 



^/s{l - s) 

1 — s 

—s 

- s] 



\ l-s 

-^s{l - s) 

-x/s(l - s) 

s 

/ 


( "" 

Vs{l - s) 

-\/s(l - s) 

l-s 

\ 


^ys{l-s) 

s 

l-s 

^5(1 - s) 



-VsC^ - s) 

l-s 

s 




\ l-s 

^5(1 - s) 

-\/s(l - s) 

—s 

/ 


( 

-Vs(l - s) 

^/s{l - s) 

l-s 

\ 

_l_ 

-^s{l-s) 

s 

l-s 

-V's(i - s) 


zt 

^5(1 - S) 

l-s 

s 

^5(1 - s) 



\ l-s 

- s) 

^5(1 - s) 

—s 

/ 


( 

-Vs{l - s) 

-Vsil - s) 

l-s 

\ 

_l_ 

-^s{l-s) 

—s 

l-s 




-^5(1 - S) 

l-s 

— s 




\ l-s 

^5(1 - S) 

\/^(l - s) 

s 

/ 


Risis) = zb 


Now according to (1-5), we identify each column of these matrices with a 2 x 2 matrix bk', it turns 
out that we can always relate the resulting b with S by the conjugacy relation (1.16) for some P. For 

61 62 


instance, Rio{s) gives the ^-symmetric basis b = 
hi = 
b3 = 


with 


S \/s{l - s) 

y/'SCl - s) 1 - ■s 

/ y/sCl - S) I- S 


63 64 

( y/sll - s) 


-\/S 


s)J ’ 


62 = 


64 = 


—s 


V i-s -v's(i - s) 

1 - s - s) 

-v'sll - s) S 


16 



and moreover, b = ^PgbPg, with 


y/s y/1- s 
VI - S -y/s 



If we do no more allow for the action of Gq, we can look for a h-symmetric basis associated with an 
isometry R G 0_(4), its product with an element of Go being again a rotation. This readily gives 
(1.16) and Proposition 1.4 is proved. 
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